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Abstract 

We derive a general formula for the correlation function of two identical particles 
with the inclusion of multiple elastic scatterings in the medium in which the two 
particles are produced. This formula involves the propagator of the particle in 
the medium. As illustration of the effect we apply the formula to the special case 
where the scatterers are static, localized 2-body potentials. In this illustration both 
i?g ide and i?o Ut are increased by an amount proportional to the square of the spatial 
density of scatterers and to the differential cross section. Specific numbers are used 
to show the expected magnitude of the rescattering effect on kaon interferometry. 


1 Introduction 

The method of two-particle intensity interferometry was developed by Hanbury-Brown 
and Twiss (HBT) in the 1950’s as a means of determining the dimension of distant 
astronomical objects [1]. It was first applied in subatomic physics by Goldhaber et al. [2] 
to study the angular correlations between identical pions produced in pp annihilations. It 
was gradually realized that the correlation of identical particles emitted by highly excited 
nuclei are sensitive not only to the geometry of the system, but also to its lifetime. 
Nowadays HBT interferometry is the primary method to probe the space-time structure 
and the dynamic evolution of the hot matter created in relativistic heavy ion collisions, 
where a quark gluon plasma is expected to be formed [3]. A lot of work has been done 
both in the theoretical [4]-[6] and experimental [7]-[9] sides (for reviews see [10]-[12]). 

HBT measurements are often performed with charged particle pairs, mostly pions. 
They suffer the long-range Coulomb interaction between themselves and with the remain¬ 
ing source which is created by the primary protons in the colliding nuclei. As a result of 
that, corrections have to be made to subtract the Coulomb final state interactions before 
the measured correlations can give useful information about the emitting source [13]-[16]. 
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Additionally, rescattering after production is another effect that may distort the cor¬ 
relation function obtained from experiments. Although rescattering effects are negligible 
for the applications in astronomy and elementary particle collisions (ep, pp and pp colli¬ 
sions), it may be significant in heavy ion collisions. In heavy ion collisions, the number 
of participating nucleons is large (around 400 for central collisions at Relativistic Heavy 
Ion Collider (RHIC)) and the produced particle multiplicity is in the thousands. The size 
of the fireball is so big that the produced particles may be scattered a few times in the 
medium before flying off to the detectors. The usual interpretation of HBT interferometry 
is that it measures the distribution of the last scattering (freeze-out distribution). What if 
the last scattering is very soft? How could that possibly affect the two-particle correlation 
function for identical particles? This is one of the questions we address. 

This paper is organized as follows. In Sec. 2 we derive a general non-relativistic 
formula for the correlation function with the inclusion of rescattering in the medium. In 
Sec. 3 the generalization to relativistically moving fermions and bosons is worked out. 
In all cases the resulting correlation function involves the propagator of the particles 
moving through the medium to the detector. In Sec. 4 we perform a calculation in one 
dimension with the result that rescattering does not change the correlation function. We 
also perform a calculation in three dimensions for a static distribution of scattering centers 
and weak 2-body potentials. The rescattering effects on the correlation function and on 
HBT radius parameters are given. The analytical formulas for R out and R s id e are derived 
since their combination will give us information on the spatial and temporal distribution 
of the produced matter. Both R^ ut and Rg ide are increased by an amount proportional to 
the square of the spatial density of scatters and to the forward differential cross section. 
We illustrate these results from by plotting the R 0 ut/-R S ide ratio using the experimental 
parameters from RHIC for kaon interferometry. In Sec. 5 we summarize our results and 
discuss future applications. 

2 Nonrelativistic Analysis 

In this section we first develop the basic formalism for the two identical particle correlation 
function. The particles are produced with a localized wave-packet, are subsequently 
allowed to interact with the surrounding medium, and are finally observed as momentum 
eigenstates by a detector. The limits of chaotic and coherent production are treated 
explicitly. 

2.1 Basic formalism 

We first study the detection of a single particle of four-momentum k = (/c°,k) at the 
space-time point x' = (t',x'). We assume that this particle was produced as a localized 
wave-packet described by , 0prod(^)- See Fig. 1. The probability amplitude for this to 
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Figure 1: A particle is emitted at a typical source point x of an extended source, and is 
detected with momentum k at the detection point x'. 

occur is the overlap of the original wavefunction, propagated through time from t to t ', 
with the observed wavefunction 

J d 3 x<i 3 xVo bs (k; x')G(x.', t x, t)^ piod {x ). (1) 

Here G(x', t'\ x, t) is the full nonrelativistic propagator, including all interactions with the 
surrounding medium, and ^ 0 b s (k;a; / ) is the wavefunction of the particle when it is ob¬ 
served. We assume it to be a momentum eigenstate such that x') = 

(The absolute normalization drops out of the correlation functions.) The produced wave- 
function is a superposition of plane waves centered at xp at time tp with phase x( x p)- 
That is, 

= /703* , (P;^) ei * ( " ,e "‘ |E ' < ' _ ' P, " P(I " Xrtl . (2) 

where <^(p; Xp) is the momentum space wave function. Without loss of generality it can be 
taken to be real and non-negative. The phase x( x p) can be chaotic, coherent, or partially 
coherent depending on the nature of the source. We can write the probability amplitude 
as 

v4(k; x P ; t ) = e ixM e iEkt ' j xp)e i ^ tp - p ^ p) e~ iEpt G( k, t p, t) (3) 

where 

G(k, t'; p, t) = I d 3 xd 3 x'G(x', t'- x, f)e ,:p x e- ,:k x '. (4) 
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Figure 2: A particle of momentum k\ is detected at x\ and another identical particle with 
momentum k 2 is detected at x' 2 . They are emitted from the source points X\ and x 2 . 


The single-particle momentum distribution P(k), which is the probability for a parti¬ 
cle of momentum k to be produced from the extended source x and arrive at the detection 
point x', is the absolute square of the total probability amplitude, 


P(k) = 


iEw r d3 p 


e ix(x P ) e iE k t: 

Xp 

Y,^ Xp) K{k; Xp ;lf,t) 

Xp 


/ (2^^ (p; xp)e i{Eptp ~ P ' Xp) e~ iEpt G ( k , t) 


( 5 ) 


Here we define a function 

K( k; xp] t', t ) = e iEkt ' I k, t'; p, t) (6) 

as a folding of the propagator with the momentum-space wave function of the produced 
wave-packet. 

Now we turn to the detection of two identical particles. Consider a particle with four- 
momentum k\ detected at x\ and another with four-momentum k 2 detected at x' 2 , as shown 
in Fig. 2. The probability amplitude for this occurrence is the product of A(k\, Xip]t',t) 
and v4(k 2 ; x 2 p] t 1 , t). Because of the indistinguishability of the particles, the probability 
amplitude must be symmetric with respect to the interchange of two particles if they 
are bosons or if they are fermions in a spin-singlet state, and antisymmetric if they are 
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2.2 Chaotic source 

A chaotic source is one whose phase x( x p) associated with the source point xp is a 
random function of the source point coordinate. To take advantage of the randomness of 
the phases, we expand the right-hand side of Eq. (5) into terms independent of x( x p) 
and terms containing x( x p)- We thus obtain 

P(k) = ^2 |/l(k; xp] t', f)| 2 + £ e ! ^)e-^)/f(k;rp;t'i)r(k;|/p;t'd). (9) 

x p xp^yp 

The second term in the above equation gives a zero contribution because the large number 
of terms with slowly varying magnitudes but rapidly fluctuating random phases cancel 
out one another in the sum. Therefore, the single-particle distribution becomes 

-P(k) = \K(k;x P ;t',t) = j d A x P p(x P ) K(k; x P ; t', t) (10) 

Xp J 

for a continuous source, where p(xp) is the density of the source points per unit space-time 
volume at the point xp. 

Similarly, we can rewrite the two-particle distribution function P(ki,k 2 ) using the 
random nature of the chaotic phases. We obtain 

P(ki,k 2 ) = \ ^2 K(k 1 ;x 1 p)K(k 2 ; x 2p) ± K(k 1 ]X2p)K(k 2 -, x ip) 

^ XlP,X2P 

+ I ^ Li(x(xiP)+x(x 2P )-x(yi p)-x(y*p))[ K (k 1 -,x 1P )K(k 2 -,x 2 p) 

{xipX2p}^{yipy2p} 


5 



± K{k l] x 2P )K{k 2 -,x 1 p) K*{k l] y 1 p)K*(k 2 -,y 2 p)±K*(k 1 -,y 2 p)K*(k 2] y 1P ) 


= P(k 1 )P(k 2 )± 


K(k 1 ]x P ; t', t)K*{ k 2 ; x P ; t\ t) 


X P 


(ii) 


Here the only surviving terms are those with X\p = y\p and x 2 p = y 2 p , or x\p = y 2P and 
x 2 p = yip. The correlation function C 2 (ki,k 2 ) becomes 


C2(ki,k 2 ) — 1± 


d i x P e' q ' xr k 2 ; x P ) 


( 12 ) 


where q = k\ — k 2 is the difference of the two detected four-momenta. This involves the 
effective density function p e ff(ki, k 2 ; xp), which is 


p ef r(ki, k 2 ; xp) 


K(ki,xp-,t',t)e tkl ' Xp K*(k 2 ,xp\t',t)e lk2 ' xp 


P( kr 


P(k 2 


■P(xp) ■ (13) 


This is the essential result of our paper. It includes the effect of rescattering of the particles 
as they travel through the medium on the intensity correlation function. Although the 
derivation was nonrclativistic, the same formula will be shown to apply when the particles 
move relativistically, with the appropriate change of the propagator. 

It is useful to check this formula in the limit of freely propagating particles. The free 
nonrclativistic propagator in position space is 




m 

2n i{t' — t) 


exp 


im{x.' — x) 2 
2 [t' - 1) 


9(t' - 1) 


(14) 


and its spatial transform is 

G 0 (k, t'; p. t) = (27r) 3 <5(k — p )9{t' — t)e~ lEk( - t '~ t ' > . (15) 


This gives the standard result 


/, , v(ki,xp)v(k 2 ,x P ) 

Peff(ki, k 2 ; xp) — - =— P{Xp). 

/>n<, j/'tkpi 


(16) 


When the produced wave packet in momentum space is independent of production point 
then p eff = p. 


2.3 Coherent source 

An extended coherent source is described by a production phase x( x p) which is a well- 
behaved function of the source point coordinate xp. The phases x( x p) and x(up) at two 
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distinct source points xp and yp are related. A simple example of a coherent source is 
the one in which the production phase x( x p) is a constant. A second example is the case 
when the production phase x( x p) is a simple function of time, independent of the spatial 
coordinate. The two-particle momentum distribution P(k 1; k 2 ) for bosons can be written 
as 


P{ k 1? k 2 ) = 


E 


J(x(xip)+x(x2p)-x{yip)-x{y2p)) 


x 


xip,x2p,yip,y2p v 

K( k l5 x 1P )K( k 2 ; x 2 P )K*( k i; y 1P )K*(k 2] y 2P ) 

+K(k 1 -,x 1P )K(k 2 ] x 2P )K*(k 1 ; y 2P )K*( k 2 ; y 1P ) 
+K(k l -x 2P )K(k 2 -x 1P )K*( k l5 j/ip)Ar*(k 2 ; j/ 2 p) 

+X(k i; z 2P )iF(k 2 ; x 1P )K*{ k l5 j/ 2 p)lT(k 2 ; j/ 1P ) 


(17) 


All four terms are essentially the same because all the production points are to be summed 
over. Obviously this is exactly equal to P(ki)P(k 2 ). Hence the two-particle correlation 
function C 2 (ki,k 2 ) = 1 for a coherent source. The probability for the detection of one 
boson is independent of the probability for the detection of another boson. The sym- 
metrization and the factorization of the probability amplitude insure that there is no 
correlation for a coherent source. 


3 Relativistic Analysis 


In the relativistic domain, the generalization is quite natural. For example, spin-0 bosons 
with mass m moving under a one-body optical potential V(x) satisfy the Klein-Gordon 
equation [5] 

(d^dn + m 2 + V(x))(f>(x) = 0 . (18) 

The propagator iA v (x',x) = (O|T(0(a/)0(a;))|O) satisfies the integral equation 


A v{x'i x) = Aq{x' - x) + J d 4 yA 0 (x' - y)V(y)A v (y, x) (19) 


with A 0 being the free Feynman propagator. Naturally, Ay reduces to A 0 in the case 
V —> 0. The general solution to this integral equation involves a summation of positive 
and negative frequency solutions, that is, 4>{x) = cj)^ + \x ) + <p^\x). We consider only the 
positive frequency part which is propagated forward in time by 


(j^ + \x') = — j d 3 x 


A v (x',x) do 4> i+ \x) 


( 20 ) 


where a do b = a ( db/dt ) — ( da/dt ) b [17]. 
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In analogy with the nonrelativistic case, the probability amplitude that a boson is pro¬ 
duced at x and detected at x' is the overlap between the propagated initial wavefunction 
with the observed wavefunction 


where 


A = - j d 3 xd 3 x'<j)* ohs (k,x') A v (x',x) do </> P rod(x ) 


0 prod (a;) = [ ^ <p(p ,x P )e ix ^e- ip < x -^ 

J J 9R 


2E p (2n 


<l>ob s (k,x') = 


—ik-x' 




( 21 ) 


( 22 ) 

(23) 


The normalization factors do not matter since they will cancel out in the correlation 
function. The K function defined in Eq. (6) changes to 


K (k; xp; t', t) = — d 3 xd 3 x 


3_/ d3 P <P{P,X P ) ip-(x-xp) -ik-x' 


(2vr) 3 J4E k E p 


A v (x',x) do e~ iE ^- tp) 


(24) 

The effective density p e g and correlation function can be calculated if the explicit form of 
the propagator is given. 

It is easy to evaluate the correlation function for free relativistic bosons simply by 
using the free propagator Ao- One obtains exactly Eq. (16); this is left as an exercise for 
the reader. 

Similarly, relativistic spin-1 fermions with mass M obey the Dirac equation 


( p — M — V(x))'ijj(x) = 0 . 


(25) 


Unlike bosons, the intrinsic spin degree of freedom can distinguish the same kind of 
fermions with different spins. The spins of the observed fermions need to be specified 
explicitly. The positive frequency solution is now propagated by S' F (x',x), which is 


S' F (x', x) = S F (x' — x) — i J dt S F (x' — x)V (x)7 0 


(26) 


where S F {x' — x) is the free Feynman propagator from x to x'. The amplitude is 

A = i J d 3 xd 3 x'^l hs (k,x\s')S' F (x\x)^ vroA (x). (27) 


The observed wavefunction is 

V'obs (k,x',s') = 


1 ^e- ik - x 'u(k,s') 


(2(/r) 3 / 2 V E k 


(28) 



where u(k, s') is the spinor which is a positive-energy solution of the Dirac equation with 
momentum k M and spin s'. Moreover 




.. ±1/2 - < 2”-) 3/2 

with b(p, s ) being the momentum space wavefunction. The K function changes to 


,3 ,3 / 

a xd x 


A (k; xp; / ,/) — ^ TTTTs M hTTh 

^S/2- 7 (2^) 3 \| EkEp 

xu\k, s')S' F (x', x)^q u(p, s ). 


M 2 


b(p, s)e 


-i\s.’x!^—ip-(x—xp) 


The spin should be taken care of in the antisymmetric two-particle amplitude, 

1 

•'ttwo particles "^= [ 4 (/ji, 3 qp, S^)* 4 .(/u 2 j ^2) A(k 1, X2P, <Si )M(A)2, Xip, S2)] 


(29) 


(30) 


(31) 


The correspondence between relativistic and nonrelativistic fermions is the substitu¬ 
tion of iS' F (x',x )70 for G(x',x). For free relativistic fermions one obtains 


Grf ree fermion — Id X(l X [/' S p( X , e e 


kP X Q“* kx ' 


= i(27r) 2 h(k - p) / dq 0 


(g 0 7° - k • 7 + M)7q 
9o - El + ie 


(32) 


In the nonrclativistic limit this reduces to Eq. (15). 

Comparing with the results in Sec. 2 , we find that the difference between nonrcla¬ 
tivistic and relativistic formulas is just the change of propagators. The nonrclativistic 
propagator G(x\ x) transforms to 


for relativistic bosons and 


G(x',x) —> -A v (x',x) do 


G(x', x) —> iS' F (x', x)j 0 


(33) 


(34) 


for relativistic bosons and fermions, respectively. As expected, the amplitudes in Eq. (21) 
and (27) both reduce to Eq. (3) in the low-energy limit, as already mentioned. 


4 Rescattering from Weak Potentials 

The usual interpretation of HBT interferometry is that it measures the space-time dis¬ 
tribution of the last scattering, the so-called freeze-out distribution. However, if the 
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scattering is very soft and clastic, we should go backwards in time and take into account 
multiple scattering effects. Then the boson propagator G(x', t'; x, t) will become 


£?o(x', t'-, x, t) + Gi(x / , t'; x, t) + G 2 (x', t!\ x, t) H- 

where G n (x.', t'; x, t) is the propagator with n scatterings [17]. In 3-dimensions we have 


G' 0 (x / ,f';x,t) = 


m 

2 

im(x' — x) 2 ] 

2m{t' — t ) 

GXp 

-1 

1 

CM 

_i 


- t), 


(35) 


Gi(x',t';x,t) = —ij d 3 x s dt s G 0 (x r ,t'; x s , t s )V"(x s , t s )G 0 (x s , t s ; x, t), (36) 

G 2 (x',t';x,t) = H ) 2 / d 3 x ls dt ls d 3 x 2s dt 2s G 0 (x / ,t';x 2s ,t 2s )V{x 2 s,t 2 s) 

x Go(x 2s , t 2s ] x is , ii s )V(x ls , t ls )Go(x ls , f i. s ; x, 0 • (37) 


For sake of illustration we take the 2-body potential to be both time-independent and 
a function only of the distance between the scattering center, located at x c , and the 
scattering point, located at x s . 

V (x s , t a ) = V(x s - x c ) = J ^-^e* q (Xs " Xc) h(q). (38) 

Here u(q) is the Fourier transform of the 2-body potential. For a Gaussian potential 
u(q) = u(0)e -q r o and for a Yukawa potential u(q) = u( 0 )/(l + q 2 rQ). A weak potential is 
characterized by a small value of u(0). We assume that the separation between scattering 
centers is large compared to the range r 0 . We also assume small momentum transfer for 
which |q|r 0 <C 1, so that a Gaussian representation is adequate. 


4.1 Rescattering in 1-dimension 

As a first illustration consider the one-dimensional problem. A particle produced at 
xp = 0 and observed at x = oo has many possibilities for scattering from localized 
potentials scattered along the x-axis. The particle can propagate directly to the detector. 
Or it may scatter once from a potential localized along the positive x-axis, or it may 
back-scatter from a potential localized along the negative x-axis. The former involves 
a transmission coefficient, the latter a reflection coefficient. It may scatter twice: from 
potentials both along the negative x-axis, or from one along the negative and the other 
along the positive x-axis, or from both along the positive x-axis. When the transmission 
coefficient is much greater than the reflection coefficient (high energy and weak potentials) 
it is easy to show that 

A'/A'free = 1 — i^Txmv^S)/k — ( 2n) 2 m 2 v 2 (0)/2k 2 + • • • = exp(— i2nmv(0)/k). (39) 
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Since this is independent of the production point Xp the correlation function is unaffected 
by rescattering and p e fr = P- This is true even when back-scattering is included. 

In one dimension, when the energy is high enough and the individual scattering centers 
are localized, HBT measures the size of the production source where the wave-packets are 
produced. Elastic scattering on localized potentials does not affect this result. This 
is consistent with an analysis in [18] which considered high energy rescattering in 3- 
dimensions using the Glauber approach. 


4.2 Rescattering in 3-dimensions 

Next consider rescattering in three dimensions. We assume for simplicity of illustration 
that the distribution of the scattering centers is Gaussian, 


Pscat (X c ) 



(40) 


Here R c and N c are the radius and the total number of the scattering centers, respectively. 
We also assume that the source function is Gaussian and parametrized as 


p(x) = 


4n 2 RpT 

After averaging over all scattering centers, we get 

Ao(k; x P ] t', t) = e lk ' xp (p( k), 
Ad(k; x P ; t', t) = 

AT 2 (k; x P ] t', t) = 


e -x 2 /2 r 2 p e -P/2r 2 


• A / 

i —expl - 

k L 

-^Cpj2t 2 c ) 

A ' / 
2|k| 2eXP t 

-Xp±Ac} 


(41) 

(42) 

(43) 

(44) 


Here we define l 2 = R 2 + 2rQ and A = mu(0)./V c /27r/ 2 . The x PP is the projection of x P 
onto the plane which is perpendicular to the momentum k. We have also assumed that 
the momentum transfer is small compared to the detected momentum k, in accordance 
with the assumption that the 2 -body potential is weak. As can be seen from above, 


K( k; xp] t ', t) 
A' 0 (k; x P \t',t) 


exp 




(45) 


The one-particle momentum distribution P(k) is unchanged by the multiple scattering. 
Therefore, the correlation function takes the form 


C2(ki,k 2 ) — 1+ 


d 4 xp(x)e iq ' x ex p) iA 


e -L 2 


/2 i 2 c 


e-xii/2 il- - 12 


(46) 
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For the two-particle momentum distribution P(ki, k 2 ), we are interested in two special 
cases: ki || k 2 and |ki| = |k 2 |, because they will provide information about P out and -R S ide, 
respectively. Dehne pair momenta k = (ki + k 2 )/2 and q = ki — k 2 . If ki is parallel to 
k 2 , we have x_|_i = xj_ 2 and k || q. Then the correlation function becomes 


CMk^ka) 


1 + e -{^WR 2 P ) 


A 2 R%l 2 c / 1 

(p c + R 2 P y(p c + 2 ^) 11 ^ 



(47) 


For not too small total momentum |k| the source size can be determined from the curvature 
of the correlation function at |q| =0 [6]. We get 


E> 2 
-^out 


2 k 2 r 2 A 2 R%l 2 c 1 

P + m 2 + (l 2 c + Refill + 2Rp) |k| 4 ' 


(48) 


If ki and k 2 have the same magnitude, k is perpendicular to q. After the integration in 
Eq. (46), we obtain 


R 


2 

side 


2 A~Rp 1 

p {ll + 2R 2 P ) 2 |k| 4 • 


(49) 


Rescattering increases both effective radii compared to the primordial source radius, as 
should be expected, albeit by differing amounts. 


4.3 Numerical example: kaons 

As a numerical example, consider kaons scattering from slowly moving nucleons. The 
forward differential cross section is taken to be 1 mb/sr. The central density of nucleons 
is taken to be 0.15 nucleons/fm 3 and the radii are Rp = R c = 5 fm. The duration of 
emission is taken to be r = 2 fm/c. 

The outward correlation function is shown in Fig. 3. A similar plot can be drawn for 
the sideward correlation function. Since the HBT radii are inversely proportional to the 
width of the correlation function, rescattering in the medium tends to increase both HBT 
radii. The resulting value of P 0 ut/-Rside is shown in Fig. 4. We see that rescattering tends 
to reduce the -Rout/Aside ratio, especially for low-energy particles which spend a longer 
time in the medium. The ratio is less than 1 for k < 150 MeV and greater than 1 for 
k > 150 MeV. The ratio is a monotonically increasing function of k. The derivation given 
in Sec. 2 is quite general, but this numerical example does not include any collective flow 
of either the produced particles or of the scattering centers (nucleons), which would tend 
to make the individual radii decrease with increasing momentum. As a result of that, it 
is not appropriate to compare Fig. 4 with RHIC experimental data. 
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Oout(p) 



Figure 3: The outward correlation function C out (q) as a function of relative momentum q 
with the parameters specified in the text. 



k(MeV) 

Figure 4: The ratio of outwards to sidewards radii that would be measured for kaons. 
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5 Conclusions 


In this paper we have derived explicit formulas for the Hanbury-Brown Twiss intensity 
interferometry correlation function when the particles undergo rescattering in the medium 
in which they are produced. The primary result is encapsulated in Eq. (12) (for chaotic 
sources). It is essentially expressed in terms of the Fourier transform of an effective source 
density given by Eq. (13). The latter involves a folding of the propagator of the particle 
in the medium with the momentum-space wave function. For nonrelativistic particles one 
uses Eq. (6), while for relativistic bosons one uses Eq. (24) and for relativistic fermions one 
uses Eq. (30). The correspondence between the nonrelativistic and relativistic versions 
of the propagators is given by Eqs. (33) and (34). It should be emphasized that, in 
general, these propagators are time-dependent since the medium is most likely expanding 
and cooling. 

For examples, we studied the rescattering from weak potentials in the nonrelativistic 
limit. In one dimension there was no effect on the two-particle correlation function. 
In three dimensions, we found that both Rg ide an d -Rout are increased by an amount 
proportional to the square of the spatial density of scatterers and to the differential cross 
section. Specific numbers were used to show the expected magnitude of the rescattering 
effect on kaon interferometry. In particular, the -R 0 ut/-Rside ratio increases with transverse 
momentum from below unity to above unity. 

This theory now must be applied to real experiments and data on heavy ion collisions. 
The general approach is applicable from low energy fragmentation reactions to the highest 
energy collisions where quark-gluon plasma formation is expected. Complications such 
as moving scattering centers, long-lived resonances, flow, partial coherence, and so on 
must be taken into account for detailed comparison to experimental data. Such work is 
in progress. 
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